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We stndy B meson decays to two charmless baryons in the diquark model, including strong and 
electroweak penguins as well as the tree operators. It is shown that penguin operators can enhance 
R —> BaB considerably, but affect R —> B1B2 only slightly, where B(i^2) and B3 are non-strange and 
strange baryons, respectively. The 7 dependence of the decay rates due to tree-penguin interference 
is illustrated. In principle, some of the BsB modes could dominate over B1B2 for 7 > 90°, but in 
general the effect is milder than their mesonic counterparts. This is because the 0 & operator can 
only produce vector but not scalar diquarks, while the opposite is true for Oi and O4. Predictions 
from diquark model are compared to those from the sum rule calculation. The decays B ^ B3B3 
and inclusive baryonic decays are also discussed. 


I. INTRODUCTION 

B meson decays provide a unique setting for baryon 
pair production, since it is impossible for the D system. 
Once observed, these decays could shed light on our un¬ 
derstanding of baryon production, and may offer further 
probes Q of underlying weak decay dynamics such as 
CP violating phases. 

Many rare mesonic B decays have been observed at 
the 10“® level in recent years, heralding the start of the 
“B Factory” era. However, rare baryonic decays have 
yet to be discovered. The most recent published lim¬ 
its come from the CLEO collaboration 0. Based on 
5.8 million BB events, CLEO finds B Ap, ApTr~ and 
pp < 0.26, 1.3, and 0.7 x 10“®, respectively. There was 
some 2.8(7 excess in the B^ pp channel, but it was 
insufficient to claim discovery. The B factories, i.e. the 
Belle and BaBar Collaborations, have now each accumu¬ 
lated an order of magnitude more data. A preliminary 
result from Belle has pushed the B^ pp limit down 
to the 10“® level. This rules out the CLEO hint, and 
puts two-body baryonic modes in strong contrast to the 
corresponding mesonic modes. However, though still elu¬ 
sive, it is quite possible that charmless baryonic modes 
are just around the corner. 

Theoretical work on rare baryonic decays is sparse. 
Most of them were stimulated by the surprising (and 
false 0) 1987 results § of Br(H“ ^ ppTr~) = (3.7 ± 
1.3 ± 1.4) X 10“'* and Br(B° ^ ppTr+Tr") = (6.0 ± 1.3 ± 
1.4) X 10“^ from the ARGUS Collaboration. Pole mod¬ 
els and the sum rule approach Q have been pro¬ 
posed for calculating the two body decay widths, while 
Ball and Dosch pursued the diquark model appoach. 
All of these works are now a decade old. Not until very 
recently, sensing that experiment is about to move for¬ 
ward, did theorists start to pay attention again. Hou 
and Soni [|| pointed out the need for reduced energy 
release on the baryon side, e.g. charmless baryonic B 
decays may be more prominent in association with rj' or 


7. Chua, Hou and Tsai studied B —> D*~NN and 
B^ —> p~pn, Tr~pn |11| using a factorization approach of 
current produced baryons. 

Eor two body baryonic decays, the sum rule and di¬ 
quark model approaches are the most relevant. The 
sum rule calculation 0 predicts that the branching ra¬ 
tio of H ^ B1B2 is typically of 10“® order while 
Br(H ^ BaB) « (0.3 - 1.0) x IQ-^ and Br(H ^ 
BcB(c)) « 0(10“^). Here B(i 2 )j Bg and Be (fenote 
non-strange, strange and charmed baryons, respectively. 
The CLEO limit of Br{B° pA) < 0.26 x lO-'^ 
1^ is already at odds with the sum rule prediction of 
Br(H+ — !■ pA) ~ 1.0 X 10“®. Eurthermore the predicted 
Br(H° ^ pp) ^ 1.0 X 10“® is an order of magnitude 
smaller than B~ ^ pA. This is in contrast with the di¬ 
quark model 1^ which typically gives a larger rate for 
B —> B1B2 over BgB, although it can predict really 
only relative rates. For example, the B^ —> pp mode 
has larger rate than all the B ^ BsB decays. Un¬ 
fortunately, Ref. 1^ did not include penguin operators. 
Judging from the role played by penguins in mesonic 
B decays like B tttt, Ktt, an enhancement by pen¬ 
guins to Br(i? ^ BgB) is to be expected. Furthermore, 
with penguins ignored, Br(H“ —> Ap) = 0 in the di¬ 
quark model, preventing us from comparing with Belle 
and CLEO search directly. 

The purpose of this paper is to complete the diquark 
model treatment of two body charmless baryonic B de¬ 
cays by including penguin diagrams, and to assess the 
importance of these penguin effects. Whether the ap¬ 
proach of using diquarks to describe baryon formation 
is correct or not is still an open question. Even if one 
assumes the idea is reasonable, to calculate the relative 
baryonic decay rates, we still need to make many other 
dynamical assumptions, which introduce further uncer¬ 
tainties. Our goal in this paper is simply to clarify the 
actual predictions of the diquark model by expanding on 
the work of Ref. Q. The experimental measurements in 
the future will decide whether diquark model or sum rule 
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approach is more relevant in the description of baryonic 
decays, or how they might be improved upon. Thus, we 
do not attempt at improving the diquark model towards 
absolute rate calculations. 

We find that penguin operators indeed could enhance 
B BgB decay rate by a factor of ^ 5 for 7 = 90°. 
Due to tree-penguin interference, the decay widths now 
depend on the unitarity phase angle 7 (= argt^*^, in the 
convention of PDG ). The enhancement in baryonic 
B decays is milder than in the mesonic decays because the 
operator Oq can not generate scalar diquarks. Penguins 
also affect non-strange decays, B —> B 1 B 2 , introducing 
also a 7 dependence, but the effect here is small. In 
general, S ^ B^B is still smaller than B B 1 B 2 , and 
B° —> pp typically has th e larg est rate. But for large 
7 > 90°, the B S+p, E+A++ rate could become larger 
than B^ pp. We find that the pattern of decay widths 
calculated using diquark model is quite different from the 
sum rule results. More experimental data should shed 
light on the two models. 

This paper is organized as follows. We first review the 
diquark model and baryonic B decays. The connection 
between penguin and diquark operators is discussed. In 
Sec. Ill, we study inclusive baryonic decays, with two 
body exclusive decays discussed in Sec. IV, in both cases 
including the effect of penguins. The conclusion is given 
in the last section. 


II. DIQUARK MODEL AND PENGUIN 
OPERATORS 


It is well known that the strong force between two 
quarks in a color-antitriplet combination is attractive, 
hence it has been speculated for a long time that they 
will form a bound or correlated state, called the diquark. 
The flavor antisymmetric combinations form scalar di¬ 
quarks while flavor symmetric combinations form vector 
diquarks. The diquark picture is useful in the description 
of baryons. The spin 1/2 octet and spin 3/2 decuplet 
baryons can be understood as bound states of a quark 
and a scalar or vector diquark, respectively jra. 


Diquarks can be generated in weak decays 114 - ^] . The 
tree level weak decay effective Hamiltonian is 


l—U,C 

-l-h.c.. 




Ci{p)0\{p) + C2{p)Ol{p) 


( 1 ) 


where 


Oi = {qpUa)v-A{Uabi3)v-A, 

O 2 = {qu)v-A{ub)v-A, 

and likewise for 0°2- The quark q could be d or s. 
The current-current operators are defined as {qiq2)v-A — 
7/^(1 — 75 ) 92 • The operators could be rewritten, after 


a Fierz transformation, in terms of scalar diquark field 
operators. For example: 

[ciOi -b C2O2] = -(c 2 - Cl) {ub)Lkiud)lk 

-b color sextet current, ( 2 ) 


with the scalar diquark field operator defined as 

{qQ)Lk — ^klmidl (1 75)Qm); (3) 

where k,l,m are color indices. Here q^ = q^C. The 
scalar diquark field operator can create a scalar diquark 
from vacuum with the strength g^Q, usually called “di¬ 
quark decay constant”: 


(0|efei„(9f 75 Qm)\{qQ)V') = )J^SiigqQ. (4) 

We see from Eq. (^) that a b quark can decay into a 
scalar diquark plus an antiquark. Since the baryons are 
bound states of a diquark and a quark, it is natural to 
expect that in decays to baryonic final states, the diquark 
operators will dominate and the sextet current operators 
can be ignored. Following this reasoning. Ref. gives a 
picture of two body baryonic B decays. The antiquark 
produced in the decay combines with the spectator an¬ 
tiquark to form a scalar or vector anti-diquark. As the 
diquark and anti-diquark fly apart, they pull a quark- 
antiquark pair out from the vacuum, resulting finally in 
a baryon-antibaryon pair. Of course, this may not be 
the only mechanism, but it is assumed to be the dom¬ 
inant one in the diquark model [^. It is interesting to 
note that O 1.2 can only generate scalar diquarks, hence 
B decaying to a decuplet baryon plus either an octet or 
decuplet antibaryon, H ^ B B , are predicted to have 
small rates 1^. These decays can only arise from the 
penguin operators as will be discussed below. 

In addition to the tree level effective Hamiltonian, it is 
well known that penguin diagrams are important for the 
charmless decays of B mesons. In mesonic decays like 
B —> Kn and tttt, penguin diagrams are crucial in the 
calculation of decay rates and CP asymmetries. Their 
effects can be described by the effective penguin opera¬ 
tors O 3 through Oio, 


n 


penguin 


-^|v/;Pt5^c,(/r)0,(^)|+h.c., (5) 

where 

^3(5) = ^')v-A{V+A){Qb)v-A, 

^ 4 ( 6 ) = -A{V+A){qoLh(3)v-A, 

q' 

^ 7 ( 9 ) = 2 ^ ^qXQ'Q')v-\-A(V-A){Qb)v-Aj 
q' 

08(10) = -^'^^q'{q0qa)v+MV-A){qabp)v-A, (6) 
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with Os-e, O 7-10 the QCD and electroweak penguin 
operators, respectively, and q = d, s, {qiq2)v±A = 
qi 7 ^(l± 75 )q 2 - The sum over q' runs over all quark flavors 
that exist in the effective field theory. 

The penguin operators can also be written in terms 
of diquark operators. Let us consider b ^ s penguins 
(similar discussion follow for the q = d case). The 
{V — A)x iy — A) type penguin operators O 34 and Ogpo 
are similar in form to Oi^ 2 - They can be written, af¬ 
ter a Fierz transformation, as the sum of diquark oper¬ 
ators '^y,q:{q'b)L{q's)\ plus color sextet terms. Again 
the latter will be ignored. The sum includes the opera¬ 
tor {uh) L (ws)^, which is of exactly the same form as that 
obtained from Oi, 2 , as well as {db)^ (ds)]^, which is not 
present in the tree operators. Note that the q' = s piece 
(s 6 )l (ss)^ vanishes since the scalar diquark is antisym¬ 
metric in their flavor constituents. 

On the other hand, the {V — A) x {V + A) type penguin 
operators O 5 and Oq do not give rise to scalar diquark 
operators. The reason is that a scalar diquark is a quark- 
quark correlation. The operators O 5 g contains one left- 
handed and one right-handed quark field, which can not 
form a scalar combination under the Lorentz transforma¬ 
tion. In other words, the scalar diquark content of octet 
baryons implies that the operators Og^g (and the cor¬ 
responding electroweak penguin operators 07 , 8 ) do not 
contribute to the B decays to octet baryon plus either an 

octet or decuplet antibaryon: S —s- BB . This is con¬ 
trary to a significant role that they play in the mesonic 
decays such as B ^ Kir, tttt. However, Og^g could gen¬ 
erate operators that consist of vector diquarks which Oi 
through O 4 could not produce. By Fierz transformation, 
for example: 

C5{syb)v^A{q'lt,q')v+A + Ce{Saybp)v-Aiq' f3'ltJ.qa)v+A 

= -^(C6 - eg) {q'RbL)Uq'RSLTk + color sextet, (7) 
with the vector diquark field operator defined as 

{dRQLTk = eklmiqfyil - 75)Qm). (8) 

Since decuplet baryons are bound states of a vector di¬ 
quark and a quark, Og g will produce B decays to a 
decuplet baryon plus either an octet or decuplet an- 

tibaryon: H B B . As mentioned above, these 
decays can not be generated by the tree Oi ,2 and the 
penguin O 37 operators. For example, the novel channel 

B~ (five strange quarks in final state) could 

arise from B~ (ss)*(su). 

To sum it up, the effective Hamiltonian that generates 
scalar diquarks can now be collected as 


Tfdiquark ^ 


Gr 




{ub)L {ud)\ + A 2 {ub)L (us)]^ 


V2{ 

+A3{sb)L (sd)l^ + A4:{db)L (ds)^ -I- h.c. 



FIG. 1 . The inclusive decay b —> Dq. 


with the coefficients 


Ai = v:yyc2 - Cl) - v;ytb{c4 - C3 + eg - cio), 

A2 = V*,,Vub{c2 — Cl) — V^*ytb{c4 — C3 -I- Cg — Cig), 

As = -VyVtb{c4 - C3 - -Cg -I- ^Cio)) 

A4 = —V)* V)b(c4 — C3 — -Cg + ( 9 ) 

For vector diquarks, we have 

T^diquark* - ^ lB4iq'RbRr{q'yR)A 

'' q'—u,d,s ^ 

+B2{qRbL)* {qRS +h.c.'^, (10) 

with the coefficients 

1 3 3 

Bi = --VtaVtb{c3 - Cg -h -e,'C7 - -e^/cs), 

B 2 = -^VtlVtb{c5 -ce + ^e,/C7 - ^e^/cs). ( 11 ) 


III. INCLUSIVE BARYONIC DECAYS 

Before we discuss the more difficult exclusive baryonic 
decays, the diquark picture could actually give us useful 
insight on the inclusive baryonic B decays. 

In the diquark model, one postulates that baryons are 
bound states of a (scalar or vector) diquark and an an¬ 
tiquark. It is natural to expect that baryonic B decays 
proceed dominantly via the process of b quark decaying 
into a scalar diquark plus an anti-quark. Since the subse¬ 
quent hadronization process always generate at least one 
baryon, with the other antibaryon guaranteed by baryon 
number conservation, the inclusive baryonic decay rates 
can be approximated by the rates of 5 ^ Dg (Fig. 1), 
where D denotes a scalar diquark such as {ud) or (cd). 

This is the approach mentioned by Neubert and Stech 
in Refs. and 0- Though citing the result from the 
above two papers. Ref. ||] adopted a different method to 
calculate the inclusive rates. It computes the rates of the 
B meson decaying into a diquark and an anti-diquark, i.e. 
B DD^ The two body channel B DD\ in which 
both diquarks are fast moving, implies that at least two 
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fast moving baryons will be generated. Though baryons 
always appear in pairs in baryonic decays, this could be 
too restrictive for inclusive baryonic decays, since it ig¬ 
nores the possibility of slow moving antibaryons. 

The b Dq rate turns out to be rather close to the 
observed inclusive rate. Let us take a closer look. After 
the b —> Dq decay, the diquark D, the antiquark q and 
the spectator antiquark jointly form a color singlet, just 
like the three antiquarks in an antibaryon. The fast mov¬ 
ing diquark D pulls out a quark q' from vacuum to form 
a baryon, leaving behind a slow antiquark q'. The color 
configuration of g', q and the spectator antiquark is again 
just like that of an antibaryon. Since q is moving fast 
while the other two are slow, the system breaks up into 
hadrons through fragmentation. This generates all kinds 
of possible final products, but at least one antibaryon has 
to be generated due to baryon number conservation. One 
possible scenario is for q to form a fast moving antibaryon 
by pulling the two antiquarks with it. Two body bary¬ 
onic decays are just such a case. Another scenario is that 
the fast moving q captures one quark to form a meson, 
leaving behind a slower antiquark. The final products of 
the decay then consist of a fast baryon, a fast meson and 
a slow antibaryon plus possible soft mesons. The baryon- 
antibaryon pair mass would then be far below rriB- One 
could also break two strings and capture two new an¬ 
tiquarks to form a fast antibaryon with the remaining 
quarks and antiquarks combining into mesons. The final 
products would then be one fast moving baryon, one fast 
antibaryon, plus two (or more) soft mesons. In all the 
above scenarios, one baryon and one antibaryon are gen¬ 
erated. But the second scenario clearly is not included 
in the B ^ DD' picture of Ref. |^. 

We list the branching ratios of 6 —*■ Dg decays in Ta¬ 
ble I. The transition amplitude of b Dq is assumed 
to factorize into the product of the diquark decay con¬ 
stant as defined in Eq. (Q), and the quark level amplitude 
{q\{qb)\b). For comparison, we also list the correspond¬ 
ing numbers for decay from Ref. by adding up 
appropriate diquark-antidiquark decay rates. For exam¬ 
ple, the rate of & ^ {cd)u would correspond to the sum of 
the rates of —> {cd){ud), {cd){ud)*. Unlike the b —> Dg 
case, calculating the latter not only involve diquark de¬ 
cay constants and masses, it also depends on B meson 
to diquark form factors. The diquark decay constants we 
use are [|jT|], 

gud, gus = 0.179, 0.215 GeV^, 

^ ^ 0.35 GeV^ (12) 

and the diquark masses are: 

UT-cs — 0.5, 0.7, 1.7, 2.0 GeV. (13) 

As expected, the inclusive baryonic decays are dom¬ 
inated by the two charmed modes: b {cd)u and 
b (cs)c. The combined branching ratio is about 4.6%. 
Adding in the rates of the smaller modes b {cs)u and 


TABLE I. Estimate of inclusive baryonic branching ratios. 
The line separates charmed vs. charmless final states. 



b^Dq 

B ^ DD' 1] 

b —> {cd)u 

2.2 X 10“^ 

2.0 X 10"^ 

b — > (cs)c 

2.4 X 10-2 

5.8 X 10"® 

b —> {cs)u 

1.1 X 10“® 

2.0 X 10-'‘ 

b —> {cd)c 

1.3 X 10“® 

3.0 X 10-'^ 

b (ud)u 

4.8 X 10-'" 

5.2 X 10"'’ 

b —> {us)u 

2.1 X 10“® 

3.9 X 10"^ 

b (ds)d 

2.0 X 10“® 

0 


b —*■ {cd)c gives a prediction for the total inclusive bary¬ 
onic B decay branching ratio of 

Br {B baryon -|- X) = 4.8%. (14) 

This is in reasonably good agreement with the experi¬ 
mental result 

Br (R ^ baryon + X) = 6.8 ± 0.6%. (15) 

The minor deficit is to be expected in consideration of 
the possibility of decaying into vector diquarks and other 
excited states as well as other mechanisms such as cur¬ 
rent produced baryons j^,Of Though Ref. ^ quotes 
a reasonable prediction from ]1^ , their approach of sim¬ 
ply adding up the i? —> DD^ rates would have given a 
branching ratio that is too small, ~ 0.8%. This is an 
indication that B DD'is not inclusive enough. In fact, 
our discussion shows that this is rather an estimate of 
the fraction of the baryonic events where both baryons 
are energetic. 

We make some observations before turning to exclu¬ 
sive modes. For baryonic decays, the single charm chan¬ 
nel b {cd)u has roughly the same rate as the dou¬ 
ble charm channel b —>■ (cs)c since the decay constants 
gcd and gcs are equal and both have two body phase 
space. This is different from the quark level picture for 
inclusive b decays, where b ccs is suppressed by a 
factor of 3 — 5 compared to b cud because of hav¬ 
ing two massive final quarks in three body phase space. 
In Ref. H], i.e. for B DD', the difference is even 
more dramatic: b (cs)c is more than twice b —> {cd)u 
because of i? ^ D form factors. This feature of the 
diquark model can be tested by experiment. For exam¬ 
ple, one can study the inclusive decays of B ^ + X 

and B —> -|- AT. These decays arise dominantly from 

b (cs)c and b {cd)u, with the diquark (cs) or (cd) 
picking up a d quark. The diquark model would predict 
Br(i? ^ 5° -b X) ^ Br(i? —> E° -b X) in strong contrast 
to expectation that b ccs is less than b cud by a fac¬ 
tor of three or more. Br {B ^ S9X) x Br (S)l ^ S“7r+) 
has been measured by CLEG |Q. While E)! decays into 
AcTt~ with 100% branching ratio, one would need abso¬ 
lute measurements of Br (5° —> S“7r+) to perform the 
test. 
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For mesonic final states, the charmed rates are 40 to 50 
times larger than the corresponding charmless rates. In 
baryonic decays, however, b —f {cd)u is 400 times larger 
than b —> (ud)u. Part of the reason is that the charmed 
diquark decay constant is larger: gcd ‘i-gud- On the 
other hand, while penguin operators enhance charmless 
mesonic decays, similar enhancement is much weaker in 
the baryonic modes, as we will discuss in the next sec¬ 
tion. We can calculate from Table I the total inclusive 
charmless baryonic decay branching ratio from the di¬ 
quark picture: 

Br {B charmless baryons -I- X) = 8.9 x 10“®, (16) 

which is relatively small. Although this estimate is prob¬ 
ably less reliable than Eq. (14), considering the numerous 
possible modes to be discussed in detail in the next sec¬ 
tion, the largest two body decay 5° —> pp is likely below 
10 “®, considerably smaller than charmless mesonic de¬ 
cays that are typically of order 10“®. 

In view of the small two body branching ratios, it is 
possible that three body decays could be larger. In a 
calculation analogous to that of ^ D*~ph 0 , it 
was estimated that ppn should be of order 10 ® 

0, hence considerably larger than two body modes. We 
note that the mechanism advocated in Ref. |^, that of 
current produced pn pair, is not contained in the diquark 
model discussed here. 



FIG. 2. R —> pp,nn through {ud){ud). 

uncertainties of the four steps, the diquark model can not 
be expected to predict absolute exclusive rates. But the 
model may give a reasonable estimate of ratios of decay 
rates. 

It will become clear that the penguin contributions 
usually involve the same or similar matrix elements as 
the tree contributions. As a result, the matrix elements 
calculated in Ref. can be used directly in our evalu¬ 
ation of the penguin effects. Since the purpose of this 
paper is to investigate the effect of penguin operators 
in the diquark model, we do not attemp at improving 
the calculation of amplitudes in Ref. ||], expecting most 
of our conclusions to be insensitive to details. We re¬ 
fer readers to Ref. for a discussion of the methods of 
calculating the various factors. 


IV. EXCLUSIVE DECAYS 


A. Non-strange decays B B1B2 


As described above, two body baryonic decays proceed 
via b —> through diquark operators. The diquark 

captures a quark from vacuum quark pair creation 
to form a baryon. The antiquark q pairs up with the 

spectator antiquark to form an anti-diquark D , which 
then captures the antiquark from pair creation and be¬ 
comes the antibaryon. Admittedly, this is not a simple 
process compared to meson pair formation. We have seen 
that Tidiquark generates only scalar diquarks from b decay, 
and hence octet baryons, while Tidiquark* generates only 
vector diquarks and hence decuplet baryons. Octet and 
decuplet antibaryons can result from b decay mediated 

by either T^diquark or T^diquark* ■ 

We find that most decay channels involve only one di¬ 
quark operator. We shall follow Ref. 0 which calculates 
the matrix element of the operators by a decomposition 
into four components: the diquark decay constant go, 
the B meson to antidiquark form factor (D \{qb)\B), the 
quark pair creation wavefunction, and the baryon wave- 
function (a diquark and a quark form a baryon). The 
authors adopt a pole model to calculate the form factor, 
take harmonic oscillator wavefunction in the ground state 
as the baryonic wavefunction and use nonlocal wavefunc¬ 
tion for pair creation. Since the latter consists of an 
undetermined normalization factor, together with other 


Let us start with the non-strange decays i? —> B 1 B 2 , 
such as B^ —*■ pp, nfi. Taking B^ —> pp as an example, 
the decay occurs only through the ^ {ud){ud) dia¬ 
gram, as shown in Fig. 2. The other decay ^ nn is 
obtained by replacing the uu pair by dd. Note that the 
decay through —*■ {dd){dd) is impossible due to the 

antisymmetry of the constituent quark flavors in a scalar 
diquark. The decay rate can be written as 


T{B° ^pp) = |Aip X \(j>p\ {ub)kiud)l \B)f. (17) 

The constant Ai, as defined in Eq. (& and evaluated at 
the scale p = mb, is equal to 4.1 x 10^ -1-4.5 x 10“'^. 

The rate without penguins, as cited in Ref. (H, is 


\V:,Vub{c2 - Cl) X \{pp\ iub)kiud)l |B)p 

4.5 X X \{pp\ {ub)k{ud)l\B)\'^. 


(18) 


Note that the two expressions differ only in the short 
distance coefficients and they share the same matrix el¬ 
ement. The matrix element {pp\ {ub)k{ud)l. \B) will be 
taken from Ref. |^. Some interesting observations can 
be made even without obtaining the absolute value. 

The rate now depends on the unitarity phase angle 7 
due to the interference of penguins with trees, analogous 
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FIG. 3. The 7 dependence of of B —> pp (dashed), T,'^p 
(solid) and Ap (dot-dashed) rates in arbitrary units. 


to the mesonic decays [Q. This dependence is shown 
in Fig. 3, together with the 7 dependence of strange de¬ 
cays discussed later. The penguin contribution is roughly 
one tenth of tree in amplitude, and its effect is milder 
compared to B ^ tt+tt". The reason is because, un¬ 
like B —> 7 r''" 7 r“, which receives sizable Oq contribution 
through chiral enhancement, the Oq operator does not 
contribute in the diquark picture, as discussed earlier. 

Other B —s- B 1 B 2 modes exhibit similar properties. 
Their rates are listed in Table II, where all rates are nor¬ 
malized to the pp mode. Decays involving a decuplet 
antibaryon, B —> B 1 B 2 , is also possible, as mentioned in 
Sec. II. An interesting channel to search for is the mode 
5“ —!■ np. Ref. ||] points out that this decay is impos¬ 
sible in the diquark picture. This prediction is still true 
even after the penguin contribution is taken into account. 
None of the candidate decay diagrams such as the tree 
B~ {ud){uu) and the penguin B~ {dd){ud) survive 
due to the antisymmetry of the constituent quark in a 
scalar diquark. This is very different from the sum rule 
calculation ||], in which Br(i3“ ^ np) is about as large 
as Br(i3° —> pp). Thus in the diquark model, one has the 
dynamical result of Br(i?° —s- pp) ~ Br(i?° —> nn), but 
Br(i?'^ ^ pn) and Br(i?“ ^ np) vanish, as seen from 
Table II. 

The tree operators Oi _2 and penguin operators O 34 
generate scalar diquarks and produce only the B decays 
into an octet baryon plus an octet or decuplet antibaryon 


—(*) 

B ^ BB . However, could generate vector di- 

quark and hence B —> is possible. Decays ^ 

A+p, A°h, A+A+, A^A*^ would arise from the vector di¬ 
quark operator {uptbi^)*{upt.di^)*\ while B° A^A~ and 
B~ A°A+, A^AO from (dab^)*{dndL)*^■ The ampli¬ 
tudes of these decays are proportional to Bi ~ 1.8 x 10 “"^, 
which is about one twentieth of Ai ~ 4.5 x 10-3. The 
vector diquark decay constants are roughly the same as 
scalar diquark decay constants |^|: 


Pud* = 0.216 GeV^, 

Pus* = 0.245 GeV^. (19) 


Assuming that the respective form factors are also of the 
same order as that for B^ —> pp, we expect the branching 

ratio of B > B*B*' ^ to be about 0.0025 x Br (R° —> pp). 

The small rates of B ^ B B is a testable feature of 
the diquark model. Because of further numerical uncer¬ 
tainties, this type of modes are not listed in Table 11. 


B. Strange decays B —> BsB 

In mesonic B decays, strange decays like B Ktt have 
larger rates than nonstrange decays like B tttt because 
of penguin contributions, with Br {B K~tt^) « 1.88 x 
10-3 compared to Br (i? ^ 7 r- 7 r+) « 4.7 x IQ-®. One 
may wonder if the same could happen in the baryonic 
decays between B ^ B^B, such as B° An, T,~^p, S°fi, 
vs B ^ B 1 B 2 such as ^ pp. This is indeed the case 
in the sum rule calculation, which gives Br (B^ —> pp) = 
1.6 X 10-3 while Br (B° ^ S+p) = 6 x 10-®. In the 
diquark caculation of Ref. Q , Br {B^ S+p) is only 
0.15 times Br (B'^ —> pp). However, since the penguin 
operators are not included in Ref. |^, it is of interest to 
include the penguins to find the actual prediction of the 
diquark picture. 

To include penguins, we proceed just like in the dis¬ 
cussion of non-strange decays. Taking B® ^ S+p as 
an example, only one diquark diagram. Fig. 4, through 
B° ^ {us){ud) will contribute, 

r(B° ^ E+p) = IA 2 P X |(Sp| {ub)k{us)l |B)|2, (20) 


TABLE II. Relative rates of B B1B2, normalized to 
r(pp) no penguin = 1 for diquark approach [9], and r(pp) = 1 in 
sum rule approach [ 8 ]. The angle 7 is taken to be 90°. 



This Work 

Ref. i 

Ref. i 

pp 

0.84 

1 

1 

nn 

0.84 

1 

0.3 

np 

0 

0 

0.6 

pA+ 

0.28 

0.33 

0.1 

nAP 

0.28 

0.33 


pA++ 

0.63 

0.75 

0.25 

nA+ 

0.70 

0.83 




FIG. 4. B° E+p . 
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FIG. 5. B —> Ap through {ds){ud). 

The constant A 2 as defined in Eq. (H), evaluated at the 
scale mb, is equal to 1.0 x 10“^ — 2.2 x 10“^. The 

rate without penguins is 

\V:Mc2 - ci)p X |(Ep| {ub),{us)l \B)\^ 

= |1.0 X 10-3e-*^p X |(Ep| iub)kius)l \B)\\ (21) 

Again the two expressions share the same matrix element, 
which we take from Ref. and find that 

|(Ep| {ub)kius)l |R)p ~ 3.0 X \{pP\ iub)kiud)l \B)\^. 

( 22 ) 

The expression for A 2 indicates that the contribution 
from penguins is almost twice as large as the tree contri¬ 
bution in amplitude and can not be ignored. The actual 
branching ratio of B^ —> E+p will depend on the an¬ 
gle 7 (see Fig. 3). The penguin operators do enhance 
the rate significantly and Br (i?° ^ E+p) is larger than 
Br (i?° ^ pp) for 7 > 90°. However, the effects are 
milder than in the mesonic decays. Br (5° ^ E+p) is 
at most twice Br (5° —> pp), when 7 = 180°. This is 
largely because the operators do not contribute to 
the same final state, unlike the role of chiral enhance¬ 
ment and constructive interference (between O 4 and Oq) 
in B ^ Ktt. Analogous to B ^ Ktt and B tttt, large 7 
will enhance Br (R° ^ E+p) but decrease Br {B^ pp). 
Similar results hold for the decays 5° ^ An, YPn. 

There are pure penguin contributions that were not 
given in Ref. |^. The modes B~ Ap, TPp are two 
examples. Nonzero tree contribution would require the 
diquark decay channel B~ {us){uu), which is impos¬ 
sible due to the antisymmetry of constituent. However, 
these modes can be generated through the diquark opera¬ 
tor {db)k (ds)/i,, which arises only from penguin operators, 
as shown in Fig. 5. This diagram is not calculated in Ref. 
1^, but is identical to ^ An, YPn, respectively, af¬ 
ter an isospin transformation u d. Hence the rate for 
> Ap is given by 

T{B ^ Ap) = IA 4 P X \{An\{ub)k{us)l\B)\\ (23) 

The coefficient A 4 is equal to —2.2 x 10“^. Since there is 
no tree-penguin interference, the rates for B~ —> Ap, E°p 
are independent of the angle 7 , just like B K^tt^ 


TABLE III. The relative rates of 5 ^ B^B normalized as 
in Table II. The angle 7 is taken to be 90°. 



This Work 

Ref. i 

Ref. i 

PP 

0.84 

1 

1 

E+p 

0.88 

0.15 

7.5 

E°n 

0.21 

0.037 


E“n 

0.72 

0 


E°p 

0.18 

0 

3.8 

An 

0.21 

0.037 


Ap 

0.18 

0 

< 3.8 

E+A+ 

0.57 

0.10 

7.5 

E°A 0 

0.17 

0.030 


E+A++ 

1.1 

0.2 

7.5 

E°A+ 

0.080 

0.014 


AAO 

0.086 

0.015 


AA+ 

0.023 

0.004 



The branching ratios of the modes R —> BgB are listed 
in Table HI. For comparison, the sum rule results 
are also listed, which gives a strikingly different pattern. 
B~ —> E+A++ and B^ E+A+ have close or larger 
rates than B^ E+p. The 7 dependence of their rates 
are identical to that of E+p. 

The 05,6 operators could generate B B*B^ ^ via 
vector diquarks. The decays B° E*+p, E*°n, E*+A+, 
E*0A0, E*+A++ and R" ^ E*0p, E*+A++ can arise 
from the vector diquark operator (unbL)* {urSl)*^, while 
BO E*-^ and 5“ ^ E*-n, E*"^, E*o^ 

from {dRbL)*{dRSL)*^■ The amplitude of these decays 
are proportional to B 2 ^ 9.0 x 10“^, which is about one 
half of A4. As described above, Br(R° —> Ap) is pro¬ 
portional to A 4 . Assuming that the form factors are of 

the same order, we expect the rates of R ^ B*B to 
be roughly 1/4 x Br (R° ^ Ap) ~ 0.03 x Br (R° ^ pp). 
This is still an order magnitude smaller than the typical 
R decays to an octet baryon plus an octet or decuplet an- 
tibaryon. Again, since further uncertainties are invloved, 

the R ^ B*B^ ^ modes are not listed in Table HI. 

C. B,B, 

This category includes decays like R° —> AA, S°A, 
etc. The relative rates of these modes are listed in Table 
IV. The mode R° ^ S°A is more straightforward since 
only one diquark diagram, B^ (us){ud), is involved. 
The enhancement effects from penguin is very similar to 
R0_-^ E+p^ 

R° —> AA decay can arise from two diquark diagrams 
and hence is more complicated. The tree contribution is 
through R° ^ {ud){ud) decay plus an ss pair creation 
(Fig. 6 (a)). The penguin operators will enhance the Wil¬ 
son coefficient in this diagram just like in the case of 
& pp. However, there is one more contribution from 
the penguin operator: through R° —> [sd)(sd) with uu 
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FIG. 6 . AA through (a) {ud){ud) and (b) {sd){sd). 


pair creation (Fig. 6 (b)). We shall argue that the penguin 
contribution is smaller than the tree. 

The short distance coefficient for penguins, ^3 = 3.5 x 
10 “"^, is smaller than one tenth of the tree V*jVub{c 2 — 
Cl) = 4.5 X 10“^. The matrix element for Fig. 6 (b) is 
basically the same as in Fig. 6 (a), except replacing ss 
pair creation by uu. To estimate the effect of this re¬ 
placement, we can compare relative rates of {B^ —> pp) 
vs AA calculated in Ref. |^, which is about 2.6 : 1. 

Since the only difference between them at tree level is just 
in the pair creation, fixes the relative weight of ss pair 
creation vs uu. This, we expect the penguin contribution 
overall to be one fifth of the tree. We therefore ignore 
the penguin contribution in our reporting of AA rates 
in Table IV. To obtain the actual numerical value and 
7 dependence in the future, one would have to evaluate 
Fig. 6 (b) . 


V. DISCUSSIONAND CONCLUSION 

In this paper, we have discussed the effects of penguin 
operators on two body baryonic B decays in the diquark 
picture. We point out that the penguin operators 03^4 
can also be transformed into operators with diquark fields 
and hence the calculation of their contributions is very 
similar to that of the tree operators. On the other hand, 
Obfi do not generate scalar diquark operators, indicating 
that their contribution to B decays into an octet baryon 
is small. As a result, penguin operators will enhance 
significantly B ^ BgB. Though the effects may not be 


TABLE IV. The relative rate oi B BsBs. Diquark rates 
are normalized so that r(pp) no penguin = 1. Sum rule rates are 
normalized so that r(pp)sumruie = 1. The angle 7 is assumed 
to be 90°. 



Rate 

Ref. H 

Ref. i 

PP 

0.84 

1 

1 

AA 

0.38 

0.39 


E°A 

0.34 

0.059 


aW° 

0.068 

0.082 


[I] 

0 

0.11 

0.02 


AE*+ 

0.068 

0.082 



0.11 

0.02 



BtB^*^) are about 0.25%, 3% 


large enough to reverse the general relative order of i? —> 
BB and 5 —> BgB as in the mesonic decays, some modes 
do have comparable rates. For example, B^ —> 
after penguin enhancement, is larger than B^ —> pp for 
7 > 90°. The B~ —*■ E+A++ mode is even of order 30% 
larger than B^ E+p. 

The penguin operators could generate B decays 
to a decuplet baryon. These channels were predicted to 
vanish in Ref. Q since the tree operators Oi ^2 only gen¬ 
erate scalar diquarks. However, {V + A) y. {V — A) type 
penguin operators 05 , 6 , 7,8 could generate vector diquarks 
after Fierz transformation, which could form decuplet 
baryons as the final product. Their rates are nevertheless 
small. We estimate the non-strange and strange decays 
Br {B ^ B*B^*^), Br(R 
of Br (i?° ^ pp), respectively. 

The diquark model calculation of the exclusive rates 
depends on the pair creation model with an undeter¬ 
mined normalization factor. Hence, absolute rates can 
not be obtained. However, as a result of duality, the in¬ 
clusive rates is independent of pair creation model. We 
estimate the inclusive rate for baryonic decays by com¬ 
puting the rate of 6 ^ Dq. This calculation relies only 
on the assumption of the diquark model and the values of 
the diquark decay constants, without further dynamical 
assumption about form factors. The total rate we get is 
very close to the experiment result, indicating that the 
diquark model is a reasonable picture for baryon produc¬ 
tion. Actually, the theoretical value is somewhat smaller, 
leaving some room for other mechanisms. 

Since the inclusive prediction relies only on the decay 
constants in the diquark model, the agreement also indi¬ 
cates that the values of pcd and pcs used are reasonable. 
Ratio of exclusive decay rates can further check the val¬ 
ues of Pud and Pus- For example, the modes B'^ ^ pp and 
5° ^ E+p have an identical B to antidiquark form fac¬ 
tor, and they differ only in the diquark decay constants 
and CKM factors. Assuming that the transition form 
factor of B to the antidiquark (ud) is not very sensitive 
to the momentum transfer, it will cancel in the ratio of 
their rates. Similar argument applies for the pair creation 
wavefunction. The ratio can be written as : 


F(i?° ^ pp) 
F(H0^E+p) 


Vcb 


9ud 

9cd 


(24) 


allowing one to in principle test diquark decay constant 
ratios. Likewise, the ratio of B^ E+A^f to B^ 
A+A^ could test pus/pcd- 

It should be clear from the previous sections that the 
diquark picture gives rather different predictions from 
the sum rule calculation. Most significantly, we note that 
H ^ B 1 B 2 is suppressed compared to B ^ B^B, B^Bc 
and BgB in sum rule treatment j|]. The reason is that, 
as the sum rule authors claim, quark pair creation is 
mainly a soft process instead of a hard one like in the 
diquark non-local pair creation model. A soft process 
would favor heavier quarks carrying larger momentum in 

















the final product to pick up soft quarks from the vac¬ 
uum. In the sum rule calculation, therefore, the ampli¬ 
tude for producing an additional quark from the vacuum 
is of order 1 in B > BcB and Be Be but suppessed in 
B —> B 1 B 2 . Such effects are much less pronounced in 
the diquark model. As a result, B —> BgB typically is 
still smaller than B —f B 1 B 2 , even after penguins are 
taken into account. 

Another feature of the diquark model is that several 
decay modes are missing due to the antisymmetry of the 
constituent quark flavor in a scalar diquark. For example, 
there is no B~ —> np while B^ —> pp and nfi have the 
same rates. The modes B~ —s- Ap, are pure penguins 

and are smaller than B^ —>■ pp. The sum rule approach 
predicts B~ —> np, Ap, S^p are of the same order as 
B*’ ^ pp. The decays arising from the penguin operator 
{ss)v±A{sb)v-A such as the novel one ^ are 

supposedly possible in sum rule calculation (though it 
is not mentioned in [^). However the above penguin 
operator {ss)v±Aisb)v-A does not have a scalar diquark 
component and thus such decays should be suppressed. 
The authors of estimate that these decays, forbidden 
by the scalar diquark model, should be suppressed by 
at least a factor of 3. The predictions emerging from the 
two pictures are, anyway, different enough to be tested in 
the near future by experimantal observation of B meson 
baryonic decays. 

The path to observation and especially understanding 
charmless baryonic B decays would be a long and winding 
one. 

Acknowledgement. This work is supported in part 
by the National Science Council of R.O.C. under Grants 
NSC-89-2112-M-002-063 and NSC-90-2112-M-003-010, 
the BCP Topical Program of NCTS, and the MOE 
CosPA Project. 


[12] D. E. Groom et al. (Particle Data Group), Eur. Phys. J. 
C 15, 1 (2000). 

[13] H. G. Dosch, M. Jamin and B. Stech, Z. Phys. C 42, 167 
(1989). 

[14] M. Neubert and B. Stech, Phys. Lett. B231, 477 (1989). 

[15] M. Neubert and B. Stech, Ann. N.Y. Acad. Sci. 578, 388 
(1989). 

[16] M. Neubert and B. Stech, Phys. Rev. D 44, 775 (1991). 

[17] B. Stech, Phys. Rev. D 36, 975 (1987). 

[18] M. Jamin and M. Neubert, Phys. Lett. B238, 387 (1990). 

[19] B. Barish et al. (CLEG Collaboration), Phys. Rev. Lett. 
79, 3599 (1997). 

[20] X.G. He, W.S. Hon and K.C. Yang, Phys. Rev. Lett. 
83, 1100 (1999); W.S. Hon and K.C. Yang, Phys. Rev. 
D 61, 0730 14 (2000); W.S. Hon , J. G. Smith and F. 
Wurthwein, hep-ex/9910014 . 


[1] W. S. Hon and A. Soni, Phys. Rev. Lett. 86, 4247 (2001). 

[2] T. E. Goan et al. (CLEO Collaboration), Phys. Rev. D 
59, 111101 (1999). 

[3] K. Abe et al. (Belle Collaboration), BELLE-CONF-0116. 

[4] C. Bebek et al. (CLEO Collaboration), Phys. Rev. Lett. 
62, 8 (1989). 

[5] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. 
B209, 119 (1988). 

[6] N. Deshpande, J. Tampetic and A. Soni, Mod. Phys. 
Lett. 3A, 749 (1988). 

[7] M. Jarfi et al., Phys. Rev. D 43, 1599 (1991). 

[8] V. L. Chernyak and 1. R. Zhitnitsky, Nucl. Phys. B345, 
137 (1990). 


[9] P. Ball and H. G. Dosch, Z. Phys. C 51, 445 (1991). 


[10] C.K. Chua, W.S.Hou and S.Y. Tsai, |hep-ph/0107110| 

[11] C.K. Chua, W.S. Hon and S.Y. Tsai, |hep-pli/010806F 


9 







